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Abstract 

The status of Lattice QCD is reviewed with respect to results that are relevant to 
Standard Model phenomenology. I argue that in a few simple cases all (or almost all) 
systematic errors from the lattice calculation are under control. Lattice QCD calculations 
of quarkonia (bb and cc) have yielded precise determinations of the strong coupling, a s , and 
the heavy quark masses. Recent work on the kaon B parameter has also led to a significant 
reduction in systematic errors, resulting in a precise determination of this hadronic matrix 
element. I discuss recent progress in other calculations of potential phenomenological 
significance, as well as problems and prospects. 



*Invited talk presented at the International Symposium Physics in Collisions '94, Tallahassee, Florida, June 
15-17, 1994. 



1 Introduction and Motivation 



Since this appears to be the first review of lattice QCD in this series, I would like to give a 
general perspective before I start the discussion of this specific review. 

By now, we have accumulated a large body of circumstantial evidence that QCD is the 
correct theory of the strong interactions. What is sorely missing is a first-principles under- 
standing of the non-perturbative effects QCD gives rise to, the most dramatic of them being 
the observed hadron spectrum. On the practical side, this lack of understanding limits the ex- 
traction of Standard Model parameters from experimental measurements. Lattice field theory 
offers a systematic first principles approach to solving QCD. Calculations in lattice QCD can 
be performed using Monte Carlo methods. 

Beyond QCD, lattice field theory has been used to study many other interesting topics, 
such as electroweak symmetry breaking, chiral gauge theories, quantum gravity, and statistical 
systems, none of which I will discuss in this review. 

I will limit the discussion in this review to calculations directly relevant to the extraction 
of Standard Model parameters. Needless to say, the field of lattice QCD is much richer than 
my limited discussion would indicate. 

Most of the work of phenomenological relevance is done in what is generally referred to 
as the "quenched" (and sometimes as the "valence") approximation. In this approximation 
(see section 2.1) gluons are not allowed to split into quark - anti-quark pairs (sea quarks). The 
obvious question from the non-expert is why such an arbitrary approximation is used at all 
in lattice QCD calculations, if the purported goal is a first-principles understanding of QCD ? 
The answer comes in two parts. Perturbative arguments suggest that most of the effect of sea 
quarks can be absorbed into the coupling constant. Of course, non-perturbative effects need 
not be so benign, and there certainly are physical quantities for which they are not. However, 
recent calculations of the light hadron spectrum (see section 4) indicate that after taking the 
continuum limit, quenched QCD can give us an acceptable phenomenology. The second part 
of the answer is calculational convenience. Calculations in the quenched approximation are 
simple enough to allow the study and control of other systematic errors that arise in lattice 
QCD calculations as detailed in section 2. They include, but are not limited to, finite lattice 
spacing, finite volume, and the extrapolation to physical light-quark masses. 

Algorithmic progress in reducing the computational expense of including sea quarks to- 
gether with the obvious increases in computing power (of which a significant part has been 
dedicated to this problem) is producing a few calculations that partially include sea quarks, 
with some phenomenological impact, as discussed in sections 3 and 5.1. 

Two limits of QCD exist where the theory becomes simple enough to be tackled with 
analytic methods. These are the chiral and the heavy-quark limits. Both limits can be studied 
in lattice QCD. In particular, the effect of sea quarks can be analyzed using techniques developed 
in the continuum: chiral perturbation theory and non-relativistic potential models. 

I would like to emphasize that lattice QCD is not a model. With all systematic errors 
under control, a lattice QCD result becomes a prediction of QCD. 

The rest of this review proceeds as follows: Section 2 gives a brief introduction to lattice 
QCD methods without going into technical details, since a number of pedagogical introductions 
and reviews already exist in the literature [1]. Attention is focused on the sources of systematic 
errors that enter lattice QCD calculations with the aim to enable the non-expert to judge the 
reliability of a lattice QCD result of interest to her/him. The following sections 3-5 discuss 
physics results by topic. Section 6 finally concludes with some remarks about future prospects. 
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2 A Lattice QCD "Primer 



Starting with the Feynman path integral formulation in Euclidean space, the discretization of 
space-time (with lattice spacing a) regulates the integral at short distances or in the ultravi- 
olet. A finite volume (of length L) is necessary for numerical techniques and also introduces 
an infrared cut-off or momentum-space discretization. The vacuum expectation of a Greens 
function, (?, which is a product of gauge and fermion fields, is defined as: 

{G) = lim lim(0) L , a , {G) L ,a = Z£ 1 a fv^VUge- 3 . (I) 

Zl,<i normalizes the expectation value. I have omitted spin and color indices for compactness. 
The gauge degrees of freedom are written as (path ordered) exponentials of the gauge field, A^. 

which makes it easy to maintain gauge invariance. The link fields, U, are SU{3) matrices. The 
(Euclidean) QCD action, 

S = S g + S f , S g = ^ jd 4 xF^F^ } S f = jd 4 xifj(x)(p+m)TP(x). (3) 

is discretized, such that Eq. (3) is recovered in the the continuum (a — > 0) limit: 

S lat = S + 0(a n ) , n>l . (4) 

I will not go into the explicit formulations of Si at here, but instead refer the reader to Ref. [1]. 
The most common form for the gauge action is Wilson's [2], written in terms of plaquettes - 
products of U fields around the smallest closed loop on a lattice. For fermions the situation is 
more complicated. The discretization of 

M = p+m , (5) 

is a sparse, finite dimensional matrix. Two different approaches are in use. In Wilson's for- 
mulation [3] chiral invariance is explicitly broken, but restored in the continuum limit. The 
pay-off is a solution of the so-called fermion doubling problem. Staggered fermions [4] keep a 
U(l) chiral symmetry at the expense of dealing with 4 degenerate flavors of fermions. 

Eq. (I) emphasizes that QCD is a limit of lattice QCD. However, in numerical calculations 
these limits cannot be taken explicitly, only by extrapolation. This is feasible, because theoret- 
ical guidance for both limits is available. The zero-lattice-spacing limit is guided by asymptotic 
freedom, since the lattice spacing is related to the gauge coupling by the renormalization group. 
Quantum field theories in large but finite volumes have also been analyzed theoretically [5]. 

In a numerical calculation the limits are taken by considering a series of lattices, as 
illustrated in Figure 1. While keeping the physical volume (or L) fixed, the lattice spacing 
is successively reduced; then, keeping the lattice spacing fixed the volume is increased. The 
calculation is in the continuum (infinite volume) limit once the hadron spectrum or matrix 
elements of interest become independent of the lattice spacing (volume). 

In practice, however, limitations in computational resources do not permit the ideal lattice 
QCD calculation just described. In particular, the computational cost of reducing the lattice 
spacing naively scales like (L/a) 4 . (The computational cost is really higher, because of numerical 
problems at smaller lattice spacings.) Eq. (4) illustrates an alternative. By improving the 
discretization errors in the lattice action (and operators), the continuum limit can be reached 
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Figure 1: Illustration of the continuum and infinite- volume limits. 

at coarser lattice spacings than before. Simulations with improved actions can come at only a 
slightly higher computational price. The ideas underlying improvement were developed some 
time ago [6, 7, 8], and have since been revitalized [9, 10, 11, 12]. 

If the quark mass is large compared to the typical QCD scale, Aqcd, effective theories 
are most adequate in describing the physics [13]. In that case, the lattice spacing cannot be 
taken to zero. Lattice-Spacing errors can, however, be systematically reduced by improvement 
[14]- 

Typical lattice QCD calculations have lattice spacings in the range 0.05 fm < a < 0.2fm, 
and finite volumes with spatial lengths in the range 1.0 fm < L < 3.0 fm. 

Now that the problem is (more or less) set up, I will briefly go through the organization 
of a typical lattice QCD calculation in the next 3 subsections. 

2.1 Gauge Configurations 

The Euclidean path integral has a formal analogy to statistical mechanics. This analogy has 
been exploited in adapting tools, originally developed for statistical mechanics systems, to 
lattice QCD. 

Because the fermion fields anti-commute (they belong to a Grassman Algebra), their 
integration can be done analytically, leaving only the integral over the gauge fields to be done 
numerically. For example, the normalization Z, becomes 

Z L:a = jvU det M[U] e~ Sa[u] . (6) 

The only viable numerical technique for the evaluation of such a high-dimensional integral is the 
Monte Carlo method. New gauge fields are sampled with probability DU det M [U] e~ Sg ^ using 
importance sampling. A gauge configuration is a lattice ensemble of such gauge fields. Approx- 
imating the integral over gauge fields in Eq. (1) with a finite sum over the gauge configurations 
introduces a statistical error. 
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Rather than going into details about algorithms, let me just note that the presence of 
the determinant, det M , complicates the process of choosing new gauge configurations. The 
quenched approximation [15] sets det M = 1 in Eq. (6), which reduces the computational effort 
by several orders of magnitude at the expense of neglecting vacuum polarization. In pertur- 
bation theory this can be compensated by a shift in the coupling constant. The interaction of 
valence quarks with gluons is, however, still treated correctly. 

2.2 Quark Propagators 

Calculating a quark propagator 

(^(x)^(y)) Li0 = Zllj VUAetM[U]M^[U]e- s ^ (7 ) 

reduces to the computational problem of calculating the inverse of the matrix M. M is singular 
for massless quarks, and the inversion becomes difficult (if not impossible) for small quark 
masses, depending on the volume, lattice spacing and sophistication of the algorithm used. 
In practice this has restricted calculations involving light quarks to masses m g >m s /3, making 
extrapolations to the physical m u ^ quark masses necessary. In some cases, this extrapolation 
can be guided by chiral perturbation theory. However, if the quark-mass dependence is not 
known a priori, control of this systematic error is possible only if the observed mass dependence 
is weak. 

2.3 Hadron Propagators 

Hadron masses are usually extracted [15] from 2 - point functions of the form: 

(GAx^y)) = (xAx)xt(y)) , (8) 

where the interpolating field, x( z )? nas the quantum numbers of the hadron in question. A 
simple choice for, say, a pion field is Xir( x ) = 4 , ( x )75 t l , ( x )- It 1S easily shown that the hadron 
propagator in Eq. (8) is a product of quark propagators in the background of the gauge fields: 

(G*(x,y)) = ZElJvUdetM[U]M y - x 1 [U] 7 5 M- y 1 [U] 7 5 e- s ^ (9 ) 

Projecting out a specific 3 - momentum, p, it can be shown that the hadron 2 - point function 
is a sum over all states (labeled n) with the same quantum numbers as, in this case, the pion: 

G„(p,t) = G^O)) 

X 

= T, c y Ent , c» = <oix» • (io) 

n 

The lowest exponential, E\ = E V} dominates at large Euclidean times, allowing the extraction 
of the pion mass from the exponential fall-off of G v (p = 0,t) at large times. 

The coefficients, C n , describe the overlap of the operators \ with the different states n. 
The presence of excited states in Eq. (10) introduces a systematic error into the hadron mass 
determination, which has to be controlled. A smart choice of \ minimizes the contributions 
of unwanted excited states to the sum in eq. (10), and thus allows for better control of this 
systematic error and the extraction of the ground-state mass at small time slices. For some 
mesons, the statistical errors (from the Monte Carlo integration) grow with Euclidean time, 
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which makes the ground-state extraction difficult. Examples include mesons with finite mo- 
mentum, orbitally excited states, and mesons with one infinitely heavy (static) quark. In such 
cases, it becomes necessary to carefully project out excited state contributions by constructing 
a matrix of hadron operators. Such techniques can also yield radially excited state masses. 

The coefficients C n become matrix elements if one of the operators \ 1S substituted with 
a current. Taking, for example, the axial current, A^x) = if)(x)~f^{ 5 if)(x) } would leave (0| A M |7r) 
in the ground-state contribution in Eq. (10), allowing the extraction of the pion decay constant. 

Three - point functions are needed in order to calculate matrix elements that are relevant 
for K — K mixing, semi-leptonic decays, etc., in a straightforward generalization of the above 
discussion. Let us consider the transition B — > A, governed by (A\J\B): 

G AJB (p A ,pj,t A ,tj) = J2 e*^e*^ (xa(x a )J(xj)xU0)) (11) 

U^oo e _ EA{tA _ tj)e _ EBtj{QlxAlA) {AlJlB) (12) 

Efficient techniques [16] have been developed to deal with the more complicated products of 
quark propagators that appear in the corresponding path integral (analogous to Eq. (9)). The 
calculational effort is about twice of that necessary for hadron masses and decay constants per 
momentum combination p^, pj. 

2.4 Perturbation Theory 

Lattice QCD calculations use perturbation theory in several places: 

• It guides the approach to the continuum limit. 

• Short-Distance quantities can be calculated non-perturbatively and compared to their 
perturbative expansions. Ref. [17] showed that indeed, 1-loop perturbation theory de- 
scribes most quantities considered to 3 — 5%, if a renormalized coupling like a-^ (rather 
than the bare lattice coupling) is used and the lattice spacing a<0.2 fm. 

• Matrix elements calculated with a lattice regulator have to be matched to their continuum 
counterpart by perturbation theory. 

• Because quarks are confined into hadrons, quark masses are always scheme dependent. 
Perturbation theory is used to convert non perturbatively determined lattice quark masses 
to the perturbative continuum masses such as the pole or MS masses. Similarly, the gauge 
coupling can be determined non-perturbatively using lattice QCD and converted to the 
MS scheme at large momenta (see also section 3). 

The lattice regulator breaks Lorentz (or Euclidean) invariance, which complicates per- 
turbative calculations relative to those performed with Lorentz (or Euclidean) invariant reg- 
ulators, such as dimensional regularization. This has prompted the development of compu- 
tational techniques for higher loop perturbative calculations [18]. (Numerical) techniques for 
non-perturbative calculations of renormalization constants have also been developed [19, 20]. 
Such techniques are very promising, because every time a new action or new operators are 
considered, not only must the programs be changed but also the perturbation theory has to be 
redone. 
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3 Quarkonium Spectroscopy 



Quarkonia, mesons containing a heavy quark and anti-quark, are at present the best understood 
hadronic systems. Both, the charm and bottom quark masses are large compared to the typical 
QCD scale, Aqcd- The bb and cc bound states are therefore governed by non-relativistic dy- 
namics. While the QCD potential was not known from first principles, relatively simple guesses 
for phenomenological potentials have proven quite successful in describing the experimentally 
measured bound state spectra of quarkonia [21]. 

As has been argued by Lepage [22], quarkonia are also the easiest systems to study with 
lattice QCD, and systematic errors can be analyzed using potential models. For example, 
finite- volume errors are much easier to control for quarkonia than for light hadrons; they are 
negligible if volumes with L > 1.5 fm are used. Lattice-Spacing errors, on the other hand, can 
be larger for quarkonia and need to be considered. However, the size of lattice-spacing errors in 
a numerical simulation of quarkonia can be anticipated by calculating expectation values of the 
corresponding operators using potential model wave functions. Control over systematic errors 
in turn allows the extraction of Standard Model parameters from the quarkonia spectra. 

Two formulations have been used in calculations of these spectra. In the non-relativistic 
limit the QCD action can written as an expansion in powers of v 2 (or 1/ra), where v is the 
velocity of the heavy quark inside the boundstate [13], which I henceforth shall refer to as 
NRQCD. Lepage and collaborators [14] have adapted this formalism to the lattice regula- 
tor. Two groups have performed numerical calculations of quarkonia in this approach. In 
refs. [23, 24] the NRQCD action is used to calculate the bb and cc spectra, including terms up 
to 0(mv 4 ) and 0(a 2 ). In addition to calculations in the quenched approximation, this group 
is also using gauge configurations that include 2 flavors of sea quarks with mass m q ~ \m s to 
calculate the bb spectrum [25, 26]. The leading order NRQCD action is used in Ref. [27] for a 
calculation of the bb spectrum in the quenched approximation. 

The Fermilab group [10] developed a generalization of previous approaches, which encom- 
passes the non-relativistic limit for heavy quarks as well as Wilson's relativistic action for light 
quarks. Lattice-Spacing artifacts are analyzed for quarks with arbitrary mass. 

Ref. [28] uses this approach to calculate the bb and cc spectra in the quenched approxi- 
mation. They considered the effect of reducing lattice-spacing errors from O(a) to 0(a 2 ). The 
results for the bb and cc spectra from all groups are summarized in Figures 2 and 3. 

The agreement between the experimentally-observed spectrum and lattice QCD calcula- 
tions is impressive. As indicated in the preceding paragraphs, the lattice artifacts are different 
for all groups. Figures 2 and 3 therefore emphasize the level of control over systematic errors. 

All groups use gauge configurations generated with the Wilson action leaving 0(a 2 ) 
lattice-spacing errors in the results. Since the heavy-quark actions have been [23], or can be 
improved through (9(a 2 ), an improved gauge action [6] should also be used, effectively reducing 
lattice-spacing errors to a negligible level. 

The first results with 2 flavors of degenerate sea quarks have appeared [26, 29] with lattice- 
spacing and finite- volume errors similar to the quenched calculations, significantly reducing 
this systematic error. However, several systematic effects associated with the inclusion of sea 
quarks still have to be studied. They include the dependence on the quarkonium spectrum 
of the number of flavors of sea quarks and the sea-quark action (staggered vs. Wilson). The 
inclusion of sea quarks with realistic light-quark masses is very difficult. However, quarkonia 
are expected to depend only very mildly on the masses of the light quarks. This systematic 
error has not been included yet and should be checked numerically. 

Quarkonia were, upon their discovery, called the hydrogen atoms of particle physics. Their 
non-relativistic nature justified the use of potential models, which gave a nice, phenomenological 
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Figure 2: A comparison of lattice QCD results for the bb spectrum (statistical errors only). 
-: Experiment; □: FNAL [28]; o: NRQCD (n f = 0) [23]; •: NRQCD (n f = 2) [25]; O: 
UK(NR)QCD [27]. 



understanding of these systems. This phenomenology is at present useful to control systematic 
errors in lattice QCD calculations of bb and cc spectra. However, we are quickly moving towards 
truly first-principles calculations of quarkonia using lattice QCD, thereby testing QCD non- 
perturbatively. In this sense, quarkonia are still the hydrogen atoms of particle physics. 

Precise determinations of the Standard Model parameters a s , m&, m c , are by-products of 
this work. They are discussed in the following two subsections. 



3.1 The Strong Coupling from Quarkonia 

At present, the QCD coupling, a s , is determined from many different experiments, performed 
at energies ranging from a few to 90 GeV [30]. In most cases perturbation theory is used to 
extract a s from the experimental information. Experimental and theoretical progress over the 
last few years has made these determinations increasingly precise. However, all determinations, 
including those based on lattice QCD, rely on phenomenologically-estimated corrections and 
uncertainties from non-perturbative effects. These effects will eventually (or already do) limit 
the accuracy of the coupling constant determination. When lattice QCD is used the limiting 
uncertainty comes from the (total or partial) omission of sea quarks in numerical simulations. 
The determination of the strong coupling, a s , proceeds in three steps: 

1. The experimental input to the strong coupling determination is a mass or mass split- 
ting, from which by comparison with the corresponding lattice quantity the scale, a -1 , 
is determined in physical units. For this purpose, one should identify quantities that are 
insensitive to lattice errors. In quarkonia, spin-averaged splittings are good candidates. 
The experimentally observed IP-IS and 2S-1S splittings depend only mildly on the quark 
mass (for masses between bottom and charm). Higher-Order lattice-spacing errors for 
these splittings are estimated in Refs. [26, 28] to be small. In Ref. [28] the calculation 
was performed at 3 different lattice spacings. 
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Figure 3: A comparison of lattice QCD results for the cc spectrum (statistical errors only). 
Experiment; □: FNAL [28]; o: NRQCD (n f = 0) [24] 



2. Within the framework of lattice QCD the conversion from the bare to a renormalized 
coupling can, in principle, be made non-perturbatively. In the definition of a renormal- 
ized coupling, systematic uncertainties should be controllable, and at short distances, 
its (perturbative) relation to other conventional definitions calculable. For example, a 
renormalized coupling can be defined from the non-perturbatively computed heavy-quark 
potential [31]. In Ref. [32] a renormalized coupling is defined non-perturbatively through 
the Schrodinger functional. The authors compute the evolution of the coupling non- 
perturbatively using a finite size scaling technique, which allows them to vary the mo- 
mentum scales by an order of magnitude. The strong coupling can also be computed from 
the three-gluon vertex, suitably defined on the lattice [33]. 

An alternative is to define a renormalized coupling through short distance lattice quan- 
tities, like small Wilson loops or Creutz ratios. For example, the coupling defined from 
the plaquette, ctp = — 31n (Tr £/p)/47r, can be expressed in terms of ay (or cxy^) by [17]: 



a P = a v (q)[l — 1.19a v (q) + O(ay)] 



(13) 



at q = 3.41/a, close to the ultraviolet cut-off. ay is related to the more commonly used 
MS coupling by 

a m (Q) = a v (e^ 6 Q)(l + -a v + ...) . (14) 
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The size of higher-order corrections associated with the above defined coupling constants 
can be tested by comparing perturbative predictions for short-distance lattice quantities 
with non-perturbative results [17]. 

Ref. [28] takes a 5% uncertainty associated with perturbation theory from this analysis. 
The comparison of the non-perturbative coupling from Ref. [32] to perturbative predic- 
tions for this coupling using Eq. (13) is consistent with this error estimate. 

In Ref. [26] the next-to-next-to-leading order corrections to Eq. (13) have been calculated 
numerically from the observed deviations (from 1-loop perturbation theory) in small Wil- 
son loops and Creutz ratios (up to size 3) at several very small lattice spacings. The 
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dominant perturbative error then comes from the conversion to the MS coupling, which 
is only known to 1-loop. 

The relation of the plaquette coupling in Eq. (13) to the MS coupling will very soon be 
known to 2-loops [34], significantly reducing the uncertainty due to the use of perturbation 
theory. 

3. Calculations that properly include all sea-quark effects do not yet exist. If we want to 
make contact with the "real world", these effects have to be estimated phenomenologically 
or extrapolated away. 

The phenomenological correction necessary to account for the sea-quark effects omitted 
in calculations of quarkonia that use the quenched approximation have been discussed in 
detail in Refs. [35, 36]; it is the dominant systematic error in this calculation. Similar 
ideas were used to correct for sea-quark effects in early calculations of quarkonia spectra 
from the heavy-quark potential calculated in quenched lattice QCD [37]. 

By demanding that, say, the spin-averaged IP-IS splitting calculated on the lattice re- 
produce the experimentally observed one (which sets the lattice spacing, a -1 , in physical 
units), the effective coupling of the quenched potential is in effect matched to the coupling 
of the effective 3 flavor potential at the typical momentum scale of the quarkonium states 
in question. The difference in the evolution of the zero flavor and 3,4 flavor couplings 
from the effective low-energy scale to the ultraviolet cut-off, where a s is determined, is 
the perturbative estimate of the correction. 

For comparison with other determinations of a s , the MS coupling can be evolved to the Z 
mass scale. An average [30] of refs. [35, 36] yields for a s from calculations in the quenched 
approximation: 



The phenomenological correction described in the previous paragraph has been tested 
from first principles in Ref. [29]. The 2-loop evolution of nj = and nj = 2 MS couplings 
- extracted from calculations of the cc spectrum using the Wilson action in the quenched 
approximation and with 2 flavors of sea quarks respectively - to the low-energy scale gives 
consistent results. After correcting the 2 flavor result to nj = 3 in the same manner as 
before and evolving a-^ to the Z mass, Ref. [29] finds 



in good agreement with the previous result in Eq. (15). The total error is now dominated 
by the rather large statistical errors and the perturbative uncertainty. 

Ref. [26] used results for a s from the bb spectrum with and 2 flavors of sea quarks to 
extrapolate the inverse coupling to the physical 3 flavor case directly at the ultra violet 
momentum, q = 3.41/a. They obtain a result consistent with the old procedure, but with 
smaller errors: 



The error is dominated by the (small) statistical errors, not the extrapolation (in nj) 
errors. The conversion to MS and evolution to the Z mass then gives: 




(15) 




(16) 







with an error now dominated by the unknown higher orders in eq. (14). 
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The claimed result in Eq. (17) (or Eq. (18)) is the most accurate determination of the 
strong coupling constant to date. In order to confirm this result, it is desirable that the bb 
and cc spectra be calculated with heavy-quark actions based on [10] with the same level 
of statistical precision and care with respect to systematic errors as was done in Ref. [26]. 
Apart from this, the systematic errors associated with the inclusion of sea quarks into 
the simulation have to be checked, as outlined in section 3. 

Phenomenological corrections are a necessary evil that enter most coupling constant deter- 
minations. In contrast, lattice QCD calculations with complete control over systematic errors 
will yield truly first-principles determinations of a s from the experimentally observed hadron 
spectrum. 

At present, determinations of a s from quarkonia using lattice QCD are comparable in 
reliability and accuracy to other determinations based on perturbative QCD from high energy 
experiments. They are therefore part of the 1994 world average for a s [30]. The phenomenolog- 
ical corrections for the most important sources of systematic errors in lattice QCD calculations 
of quarkonia are now being replaced by first-principles, which will significantly increase the 
accuracy of a s determinations from quarkonia. 

In a few years time, the world average for the strong coupling will be dominated by 
determinations of a s using lattice QCD. 

3.2 The Heavy Quark Masses 

Because of confinement, the quark masses cannot be measured directly, but have to be inferred 
from experimental measurements of hadron masses, and depend on the calculational scheme 
employed. 

In lattice QCD quark masses are determined non-perturbatively, by tuning the bare lat- 
tice quark mass {m l q ) so that, for example, the experimentally observed J/if> (or T) mass 
is reproduced by the calculation. For phenomenology useful quark masses are the perturba- 
tively defined pole and MS masses, which the bare lattice mass can be related to by (1-loop) 
perturbation theory: 

m Pole = Z pole m la t ™f (m Q ) = TU^ . (19) 

The heavy-quark pole mass can be determined alternatively from a calculation of the binding 
energy, -Ebmd- The ground-state energy need not equal the mass of a non-relativistic system. 
The binding energy can be obtained by subtracting the perturbatively calculable heavy-quark 
rest energy from the ground-state energy. The pole mass is then: 

^ = \(M^ - E hmd ) (20) 

This method is insensitive to errors in tuning the bare mass, because the binding energy depends 
only mildly on the quark mass. 

Of course, as always, all systematic errors arising from the lattice QCD calculation need 
to be under control for a phenomenologically interesting result; in particular, the systematic 
error introduced by the (partial) omission of sea quarks has to be removed. The short-distance 
corrections that introduced the dominant uncertainty to the a s determination from quarkonia 
are absent for the pole mass determination, because this effective mass does not run for momenta 
below it's mass. 
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Ref. [25] used both methods described above for a determination of the b quark pole mass 
from a lattice QCD calculation of the bb spectrum. As expected, a comparison of their results 
with zero and 2 flavors of sea quarks finds compatible results for the pole mass: 

m P° le = (5.0 ± 0.2) GeV (21) 

For the MS mass, Ref. [25] quotes m^ s (m;,) = 4.0(1) GeV. The error in both results is domi- 
nated by perturbation theory. 

A similar analysis is being performed in Ref. [38] for the charm quark mass from the 
charmonium spectrum. A preliminary result is ra^* = 1.5(2) GeV. 

The MS mass for the charm quark has also been determined from a compilation of D 
meson calculations in the quenched approximation [39], with m^ s (2GeV) = 1.47(28) GeV. 
The error includes statistical errors from the original calculations and the perturbative error. 
However sea-quark effects cannot, in this case, be estimated phenomenologically, leaving this 
systematic error uncontrolled. 



4 Light Hadron Spectrum 

The calculation of, say, the proton mass from first principles is one of the original problems 
lattice methods hoped to solve. Unfortunately, the systematic errors in the lattice calculations 
do not yet allow for first-principles calculations of the light hadron spectrum. Consequently, 
while a first-principles understanding of the (low lying) bb and cc spectra lies just around the 
corner, as I argued in the previous section, the light hadron spectrum is one of the hardest 
problems, most likely to be solved later. 




K <|> N E+E-N A 



E Q 10f„ 10f K i m 



Figure 4: The results from Ref. [41] in comparison with experiment 

Nevertheless, I would like to discuss some of the progress of the last few years [40], because 
it has been significant. The results of a systematic calculation of the light hadron spectrum 
in the quenched approximation [41] are shown in Figure 4. The inverse lattice spacing was 
determined from the mass of the p } the pion and kaon masses were used to set the quark masses. 
All the predicted hadron masses in Figure 4 are written in units of the p mass. Three sources of 
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systematic error were corrected by extrapolations in m q (to m Ut d), to zero lattice spacing, and 
to infinite volume. The comparison of the results of this calculation with the experimentally 
observed spectrum indicates that quenched QCD might not be a bad approximation to QCD. 
However, some of the assumptions used in the extrapolations still need to be checked. Lattice- 
Spacing errors can be reduced with little effort by considering an improved action instead of 
the Wilson action which was used. A similar analysis with the staggered fermion action is in 
progress [42]. 

Calculations that include sea quarks have been improved significantly over the last few 
years. However, with the exception of finite- volume effects [43], systematic errors (a — > 0, 
m q — > m u ^) have not been extrapolated away yet. On the other hand, the gauge configurations 
generated with sea quarks are useful for quarkonium physics, because lattice-spacing errors can 
be controlled and light-quark mass effects are most likely mild. 

The effects of the quenched approximation have also been studied in chiral perturba- 
tion theory [44]. This serves to guide our expectations about the limitations of the quenched 
approximation and to identify physical quantities that are insensitive to sea-quark effects. 

In principle, numerical calculations of the light hadron spectrum together with perturba- 
tion theory can lead to determinations of the light-quark masses. While this should be very 
interesting in the future, since the light-quark masses are very poorly known, the systematic 
errors at present make such determinations unreliable [45]. 



5 The CKM matrix 

The CKM matrix is a window to new physics, because it is linked to the unknown fermion mass 
generation mechanism. An overconstrained CKM matrix can give us hints about the physics 
beyond the Standard Model such as the existence of a fourth generation of quarks and leptons 
or sources of CP violation inconsistent with the one allowed phase in the unitary matrix. 
Information about CKM matrix elements can, in principle, be extracted from experimental 
measurements of exclusive weak processes, such as neutral meson mixing, semi-leptonic and 
leptonic decays, and flavor- changing radiative decays. However, these processes also receive 
QCD corrections, the long distance part of which is parametrized into a weak matrix element 
(see section 2.3, Eq. (11)). With the exception of a few special cases, these matrix elements 
cannot be calculated from first principles using analytic methods. Lattice QCD can therefore 
play an important role for this phenomenology (see also recent reviews [46]). 

5.1 K — K Mixing 

In the Standard Model K — K mixing proceeds through the exchange of two W bosons [47]: 
e ~ Im(V td )\V cb \[(f 3 (m t )r, 3 - r, iyc )\V cd \ + ri 2 y t f 2 (m t )MVu)] \m 2 K f 2 K B K , (22) 



Bk = " Lt a ' ■ (23) 



3 ( K\Oas=2 \K) 
8m 2 K f 2 K 

The fi } yi are known functions of the (charm and top) quark masses and rji are perturbative 
QCD corrections. Using three-generation unitarity, the dependence of e on the CKM angles 
can be reduced to the CP violating phase and |T4&/K&|- Bk is for a variety of reasons one of the 
simplest quantities to study in lattice QCD. Because it is the ratio of two very similar matrix 
elements, many systematic errors are expected to cancel. The chiral symmetry of staggered 
fermions makes them the method of choice for Bk [48] . 
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Most systematic effects have been studied for this quantity. While all groups see a siz- 
able lattice-spacing dependence in their results [49, 50], this error has recently been shown to 
be 0(a 2 ) on theoretical grounds [51], leading to a significant reduction of the extrapolation 
errors. The expectation that errors associated with the quenched approximation cancel has 
been checked at somewhat coarse lattice spacings with sea quarks of mass m q ~ m s /2 [50, 52]. 
The perturbative corrections to several different operators have been calculated [53, 54]. Their 
inclusion yields consistency between those operators after extrapolation to zero lattice spacing 
[50, 51]. The final result in the naive dimensional reduction scheme is [51]: 

B K (NDR, 2 GeV) = 0.616 ± 0.020 ±0.017 , (24) 

which can, of course, be converted to the more customary scale independent parameter by 
B K = a;^ 25 B K (NDR,2GeV). The first error in Eq. (24) is statistical, the second is the 
sum of all systematic errors considered. Still to be checked is the lattice-spacing dependence 
of calculations with nj ^ 0, which need not equal the nj = case. While sea quarks with 
non-degenerate masses, m q ^ ra s , still need to be incorporated into numerical simulations, 
they have been estimated [44] in chiral perturbation theory to affect Bk at the less than 5% 
level. In summary, Eq. (24) is the result of a careful analysis, with most errors under control. 
It already surpasses all other (model dependent) determinations of this quantity in accuracy. 
Once the last bit of "hand waving" with respect to sea-quark effects is checked directly, this 
result will significantly reduce the uncertainty in the Wolfenstein parameters p and rj. 

5.2 B-B Mixing 

This case is very similar to K — K mixing. The amplitude is written in the form: 

x d ~\V; d V tb \ 2 f(m t )^m 2 B t B B B . (25) 

The major uncertainty in extracting the CKM angles from the experimental measurement of 
Xd comes from the unknown decay constant Bb is expected to be close to the vacuum 
saturation value, Bb — 1 [55]. 

Heavy meson decay constants have been studied in various limits by a large number of 
groups. A lot of the simplifying features that allowed easy control over systematic errors for Bk 
are absent here. Hand waving arguments about cancellations of systematic errors can only be 
made for ratios of decay constants, /b//d, /b s //b, /d s //d, which have also been considered. 

The B meson in the static limit [56] is an example of a quantity where (for well understood 
reasons) the extraction of the ground-state signal from the correlator is difficult. This systematic 
error has been the focus of two groups [57, 58]. They use sophisticated techniques of removing 
excited states by constructing a matrix of hadron propagators over several different operators. 
Attention has also been paid to finite- volume and lattice spacing effects [57, 59, 60]. The results 
for Jb in the static limit are summarized in Figure 5. Some of the disagreement follows from 
the difficulty of extracting the ground-state signal, which has not been addressed by all groups 
in the same detail. The lattice-spacing dependence in Figure 5 is rather strong and needs 
to be better understood or reduced, before the a — > extrapolation becomes reliable. This 
dependence will be somewhat softened by the inclusion of perturbative corrections [61], which 
have been left out for the purpose of comparison. Also, using an improved action instead of 
the Wilson action for the light quarks should reduce the lattice-spacing errors. 

The decay D s — > \iv has been measured in several experiments [62], yielding determina- 
tions of fu s (assuming V cs ). Lattice calculations of this quantity can then be used to calibrate 
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Figure 5: Lattice QCD results for f^ at in comparison, o: [57]; □: [60]; O: [64]; A: [63]; V : 
[59]; x: [58]. 

the quenched approximation once all other systematic errors have been removed and the exper- 
imental results become more precise. Figure 6 shows a comparison of lattice and experimen- 
tal results. Systematic errors were considered to varying degrees and are not yet completely 
controlled (even within the quenched approximation). One group [66] has partially included 
sea- quark effects. 

For the B meson decay constant, results for mesons containing one heavy quark of order 
the charm quark (or heavier) mass have traditionally been combined with the static limit, 
and interpolated to the B meson mass. In the heavy-quark limit the combination <f> = fy/m 
becomes independent of the heavy-quark mass [68] (up to renormalization group logarithms), 
with corrections that vanish as powers of 1/ra: 

^ = ^(1 + - + ^ + ...) (26) 
ra ra z 

An example of this interpolation from Ref. [59] is shown in Figure 7; similar results were 
obtained by all groups. In particular, large deviations from the static limit as shown in Figure 7 
were seen by all groups. However, lattice-spacing errors in the heavy-quark action (and current) 
also contribute 1/ra terms. Even though they are probably small, they should be removed by 
using improved heavy-quark actions and currents [10, 14]. Improved heavy-quark methods can 
also obviate the (expensive) need to use ever smaller lattice spacings for weak matrix elements 
evaluated directly at the B meson mass. 

I conclude this section with a summary of the current state of knowledge of heavy meson 
decay constants in Table 1. I would like to emphasize again that, while attention has been paid 
to systematic errors, they are not yet under complete control in the quenched approximation. 
However, there are no technical obstacles to calculations of heavy meson decay constants in the 
continuum limit of quenched lattice QCD. While the inclusion of sea-quark effects will be more 
difficult to control, ratios of decay constants may prove to be as insensitive to these and other 
systematic errors as Bk has been shown to be. 
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Figure 6: A Comparison of lattice QCD and experimental results for fu s . The second error 
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[59]; •: ELC [64] 


O: APE [63]; 


*: UKQCD [65]; o: 


HEMCGC [66]. 




group 


f D (MeV) 


fD./fD 


f B (MeV) 


f Bs /f B 


BLS [59] 
UKQCD [65] 
PWCD [60] 
APE [63] 
ELC [64] 
HEMCGC [66] 
FNAL [57] 
(static) 
SH [67] 
(non-rel.) 


208 ± 9 ± 37 
185 ±| ±f 
170 ± 30 
218 ±9 
210 ± 15 
215 ± 9 ± 53 


1.11 ±0.02 ±0.05 
1.18 ±0.02 
1.09 ±0.02 ±0.05 
1.11 ±0.01 
1.08 ±0.02 


187 ±10 ±37 
160 ± 6±H 
180 ±50 

205 ± 40 

188 ± 23±^ 
171 ± 22±l 9 5 


1.11 ±0.02 ±0.05 
1.22±| 

1.11 ±0.02 ±0.05 

1.08 ±0.06 

1.216 ±0.041 ±0.016 



Table 1: Lattice QCD results for heavy mesons decay constants in comparison. The second 
error is the groups estimate of the systematic errors they considered. 



5.3 Semi-Leptonic Decays 

Semi-leptonic decays have traditionally been used to determine CKM matrix elements. The two 
classic examples are nuclear /3-decay for V u d and for V us . In both cases, (SU(2) or SU{3)) 
chiral symmetry together with current conservation was used to calculate the hadronic matrix 
element. Semi-leptonic decays of charm or bottom mesons have no such advantage (with the 
exception of decays like B — > D^lv, see section 5.4.1). Hence, in all but this one case exclusive 
semi-leptonic decays of D and B mesons have failed to provide accurate determinations of the 
corresponding CKM matrix elements, even though a wealth of experimental information has 
been accumulated in recent years [69]. 

The matrix element for, say, D — s- K^lv can be parametrized in terms of form factors. 
Common parametrizations are: 



(K\J,\D) = f + (q 2 



{PD +PI<) 
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Figure 7: <f> vs. 1/m from Ref. 



and 



(K*,X\J,\D) 



+ A 2 (q 2 



mj} + rriK* 
1 

I 



V(<? y mP aVDV a K* ~ M{q ){rn D + m K *)6 h 



2 2m k* 



■{PD + PK*)f,PDa ~ A(q ) 5— (j>D - PK*)fj,PD 



with momentum transfer g 2 = (p£> — Pk(*)) 2 - I n the rest frame of the D meson, zero recoil 
(Pa'(*) = 0) corresponds to maximum momentum transfer, q, 



2 

max 



(m D 



K(*), 



Maximum 



recoil is at q 



mj 



0. Knowledge of the form factors yields CKM matrix elements from 



measurements of, for example, the differential decay rate: 

dT(D -> Klv) 



dq 2 



1927r 3 m 3 D 



(29) 



or 



dT(D -»■ K*lv) 
dq 2 



IK 



r 2 \ X 2 n 2 

-{m D + m K * ) \A 1 



647r 3 m|, 



m 



D ' 



-m 



212 



-im 2 K(t) q 



q>)\> + 0(\ 3 ' 2 ) (30) 
2 It is evident from Eqs. (29) and (30) 



where A = A(ra£>, ra^(.), , 

that the rate vanishes at zero recoil, because of the kinematic factor A 1 / 2 = 2ra£>p^(.). 

It is therefore necessary to consider the form factors as functions of the recoil momentum 
(or q 2 ). However, momenta are discretized by the finite volume, and lattice artifacts increase 
as ap ~ 1. Finite- Volume errors should also be expected to change relative to the zero mo- 
mentum case. Mesons with finite momentum have similar problems with the time dependence 
of statistical errors as mesons containing a static quark (see section 5.2) or orbitally excited 
states, which makes the ground-state extraction difficult. The overall increase in statistical 
errors relative to the zero momentum case obfuscates systematic error analyses. 

All groups limit themselves to momenta with |p| < 2p m i n (p mm = ^r)- The form factors 
are calculated at the corresponding discrete values of q 2 . Smart choices of initial and final state 
momentum combinations can still cover a significant range of momentum transfers. 

The CKM angle V cs is known from three-generation unitary of the CKM matrix. The 
decay D —> K^lv can therefore be used to test lattice methods and calibrate the quenched 



16 



group 


/+(?max) 




A2/A 1 (q 2 max ) 




Exp [69] 


1.31 ±0.05 


0.66 ±0.05 


0.73 ±0.15 


2.04 ±0.27 


ELC [71] 
APE [72] 

(Ji\(alLjL) [lo\ 

LANL [74] 
BKS [70] 


1.13 ±0.31 

1.36 ±0.14 

1171 0.12 
L.U ± 014 

1.27 ±0.10 

1.64 ±0.36 ±0.36 


0.62 ±0.09 
0.79 ±0.13 

n QQ _L 0.07 

0.83 ± o 12 
0.78 ±0.04 
1.23 ±0.2 ±0.31 


0.70 ±0.40 

0.70 ±0.40 

n on -i- n on 
U.yu ± U.ZU 

0.72 ± 0.22 

0.70 ± 0.16 ± ° j° 5 


1.62 ±0.30 
1.73 ±0.32 

1 Kn -I- 0.50 

1.50 ± 20 

2.03 ±0.13 

2.15 ± 0.24 ± ° f 8 



Table 2: Lattice QCD results for D — s- K^lv in comparison with experiment. The second 
error is the groups estimate of the systematic errors they considered. 



approximation. A comparison is shown in Table 2. As pointed out in Ref. [70], another nice 
calibration mode is the semi-leptonic decay K — > 7r/z/, because precise measurements of both 
form factors, / + and / , as functions of q 2 exist. 

For historical reasons, experimentalists and (lattice) theorists alike have quoted the form 
factors at q 2 = (which is usually coupled with parametrizing the q 2 dependence by pole 
dominance), even though this kinematic point is not particularly convenient for either group. 
Eqs. (29) and (30) emphasize that experimental measurements of differential decay rates and 
lattice determinations of form factors at the same recoil momenta can be used to extract the 
CKM angles without model assumptions about the q 2 dependence. If such assumptions have 
to be made because of the experimental set-up, I suggest to quote the form factors at q 2 = ^ ax 
for comparison with theory. To emphasize this point, I have changed the standard comparison 
of form factors at q 2 = to q 2 = q^ ax in Table 2. The experimental results were simply rescaled 
by the appropriate factor ([1 — <?max/ m poie] _1 ) j f° r the lattice QCD results I took the form 
factors as quoted at ^ ax where available, otherwise they were rescaled in the same fashion as 
the experimental results. 

A study of lattice-spacing, finite- volume, and other systematic errors in Ref. [70] was very 
limited by large statistical uncertainties. While several groups have presented calculations since 
[71, 72, 73, 74], no systematic investigation of finite- volume and finite lattice-spacing effects has 
appeared yet. 

Most groups [70, 71, 72] have also considered other charm decays, like D — > Tv(p)lv, and 
D s — > ?y(</>)/z/, where possible contributions from non-spectator diagrams were neglected. 

Ref. [71, 72] extrapolated their results for charm decays to B — > ir{p)lv assuming heavy 
quark symmetry. An analysis of B meson semi-leptonic decays at realistic b quark masses has 
not been carried out yet. 

In summary, the need to consider mesons with nonzero momentum complicates lattice 
QCD calculations of semi-leptonic decays relative to calculations of decay constants. So far, 
lattice calculations of semi-leptonic charm decays have been less systematic than those of heavy 
meson decay constants. However, no serious obstacles exist, to obtaining definitive results in 
the quenched approximation at least over a limited range of recoil momenta. 

5.4 Other Weak Matrix Elements 

In the following two subsections I discuss the status of calculations of the Isgur-Wise function 
(relevant for experimental measurements of B — > D^lv) and of the exclusive radiative decay 
B -> K*-f. 

In both cases the first numerical lattice QCD results appeared only very recently. The 
lattice techniques are, of course, very similar to those established for the study of semi-leptonic 
charm decays. 
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5.4.1 The Isgur-Wise Function 



Experimental measurements of the differential decay rates for the exclusive decays B — > D*lv 
near zero recoil depend with only small corrections on the Isgur-Wise (IW) function. Using 
Heavy Quark Effective Theory (HQET) it has been shown that the corrections to the symmetry 
limit are small for this mode, yielding a model-independent determination of V c b [68, 69]. For 
present day lattice calculations this mode is useful as a calibration of the quenched approxima- 
tion. 

Two different approaches to lattice QCD calculations of the IW function [75] have been 
employed, in analogy to the B meson decay constant. The discretization of the Heavy Quark Ef- 
fective Theory allows for the direct computation of the IW function [76] from the corresponding 
matrix elements. 

The IW function can be extracted from matrix elements which are related to the above 
process by HQET. In particular, the corrections to heavy quark symmetry for the form factors 
extracted from matrix elements of the form (D\c'y l j i c\D), which were the focus of the initial lattice 
QCD calculations [77, 78], are similarly suppressed. Ref. [78] also considered the processes 
"i?"— > D^lv using heavy quark masses around charm. Figure 8 shows as an example a 
comparison of the results of Ref. [78] with experimental data from CLEO [79]. 
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Figure 8: The IW function from Ref. [73] in comparison with CLEO II data. 

The knowledge of systematic errors is similar to the case of semi-leptonic charm decays. 
No serious analysis of finite volume, finite lattice spacing and (heavy and light) quark mass 
effects exists yet. 

The error on V c b can, in principle, be reduced beyond the current knowledge, once lattice 
QCD calculations of the form factors at recoil momenta where the differential decay rate is 
measured are at least as accurate as the current theoretical estimate at zero recoil, because the 
at present necessary extrapolation of experimental results to zero recoil can then be avoided. 

5.4.2 B -> K*j 

Lattice calculations of the radiative decay, B — > K*~{ } have also recently begun [80, 81], where, 
as before, heavy quark masses of order the charm quark were used. The exclusive decay rate 
is calculated with the hadronic matrix element (K*\scr fll/ q' y bR\B) } which is again parametrized 
in terms of form factors, T\ and T 2 . The physical two-body decay takes place at momentum 
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transfer q 2 = 0. Extrapolations to mg using HQET and to q 2 = using pole dominance are 
therefore necessary. Some consistency checks are possible by interchanging the order of the 
extrapolations. 



The ratio 



„ T(B — > K*~i) r , ,o r mi,,. , , 
Rk* = -^77 r 1 = 4[m B /m 6 ] 3 [l - -JU^o)! 3 (31) 



r(6 -> s-f) 



m 



B 



has been measured [82] to Rk* = (19 ± 8) %. Lattice QCD predictions of this ratio [80, 81] are 
in the range Rk* = 6 — 9 % with an uncertainty of about 50 %. Given the exploratory nature 
of the lattice QCD calculations and the large experimental errors, this is pretty fair agreement. 

These results should be followed up by calculations directly at the B meson mass, together 
with a systematic study of lattice-spacing artifacts at large recoil. 



6 Conclusions 

Lattice QCD is becoming a reliable tool for Standard Model phenomenology, even though we 
are far from being done yet. 

With present day technology, it will be challenging to reach the level of control over 
systematic errors in calculations that include sea-quark effects that is possible in calculations 
that use the quenched approximation. This approximation will therefore be with us for some 
time. Definitive calculations of the light hadron spectrum, heavy meson decay constants and 
other simple weak matrix elements in the continuum and infinite- volume limit of quenched 
lattice QCD are possible. We will continue to see such results in the coming years. The 
comparison with the experimentally measured light hadron spectrum and weak decays (where 
the corresponding CKM matrix elements are known) can give us a calibration of quenched QCD. 
While the final judgement on whether quenched QCD can give us a useful phenomenology is 
still pending, most results in quenched QCD so far indicate that it will. 

In the meantime, progress will continue to be made in calculations including sea-quark 
effects, leading to phenomenologically interesting results and in due course to better knowledge 
of an increasing number of Standard Model parameters. The search for better lattice actions 
(with very small discretization errors at, say, lattice spacings of a ~ 0.5 fm) and better algo- 
rithms (for the inclusion of sea quarks) is a very active area of research - I invite the reader 
to take a look at the proceedings of the 1994 meeting on Lattice Field Theory [83]. This has, 
of course, the possibility of revolutionizing the field, making possible definitive results sooner 
than we anticipate. 

In some special cases, like quarkonium spectroscopy and Bk, such calculations are already 
possible with current technology and with control over systematic errors. A first-principles 
understanding of these (and hopefully other) simple physical systems should be possible within 
the next few years. 
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